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ABSTRACT 

The Pipe Nebula is a massive, nearby dark molecular cloud with a low star-formation efficiency which makes 
it a good laboratory to study the very early stages of the star formation process. The Pipe Nebula is largely 
filamentary, and appears to be threaded by a uniform magnetic field at scales of few parsecs, perpendicular to 
its main axis. The field is only locally perturbed in a few regions, such as the only active cluster forming core 
B59. The aim of this study is to investigate primordial conditions in low-mass pre-stellar cores and how they 
relate to the local magnetic field in the cloud. We used the IRAM 30-m telescope to carry out a continuum and 
molecular survey at 3 and 1 mm of early- and late-time molecules toward four selected starless cores inside the 
Pipe Nebula. We found that the dust continuum emission maps trace better the densest regions than previous 
2MASS extinction maps, while 2MASS extinction maps trace better the diffuse gas. The properties of the cores 
derived from dust emission show average radii of ~0.09 pc, densities of ~1.3xl0 5 cm -3 , and core masses of 
~2.5 M Q . Our results confirm that the Pipe Nebula starless cores studied are in a very early evolutionary stage, 
and present a very young chemistry with different properties that allow us to propose an evolutionary sequence. 
All of the cores present early-time molecular emission, with CS detections toward all the sample. Two of them, 
Cores 40 and 109, present strong late-time molecular emission. There seems to be a correlation between the 
chemical evolutionary stage of the cores and the local magnetic properties that suggests that the evolution of the 
cores is ruled by a local competition between the magnetic energy and other mechanisms, such as turbulence. 
Subject headings: ISM: individual objects: Pipe Nebula - ISM: lines and bands - ISM - stars: formation 
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1. INTRODUCTION 



The Pipe Nebula is a massive (10 4 M G : lOnishi etalJl 19991) 
filamentary (~15 pc long and ~3 pc wide) dark cloud lo- 
cated in the southern sky ~5° apart from t he Galactic Cen- 
ter. I ts short distance to the Sun (145 pc: lAlves & Franco! 
120071) places this complex in the group of nearby molecu- 
lar clouds which serve as good laboratories for star forma- 
tion surveys. Despite the large reservoir of mass, the Pipe 
Nebula molecular cloud is characterized by being apparently 
quiescent, with a ve ry low star-formation efficiency (~0.06% 
for the entire cloud. iForbrich et al.ll2CI09l) . Barnard 59 (B59), 
located at the northwestern en d of the cloud, has formed a 
small cluster of low-mass stars dBrooke et al. 2007). The low 
global star forming efficiency of the cloud contrasts with that 
of other nearby molecular clouds such as Ophiuchus or Tau- 
rus, where an important star-formation activity is observed. 
The Pipe Nebula is, hence, an excellent place to study the ini- 
tial conditions of star formation at scales of a few parsecs. 

The fi rst extensive survey toward the Pipe Nebula was 
done by lOnishi et al.l (1 1999b through single dish observa- 
tions of CO isotopologues. These authors were the first 
to suggest a clumpy distribution for the dense gas by 
detecting compact C ls O cores in the main body of the 
cloud . It was not until the last few years that several sur- 
veys dLombardi et al.1l2006MMuench et al.ll2007l:lBrooke et al.l 
120071: iRathborne et al.l 120081: lLada et al.1 120081) were car- 
ried out to explore t he physical properties of the cloud. 
lLombardi et al J (120061) use 2MASS data to construct a high 
resolution extinction map of the Pipe Nebula through which 



they identify a large number of high extinction cores with typ- 
ical masses between 0.2 and 5 M Q , Molecular line observa- 
tions reveal that they are starless cores in a very early evo- 
lutionary stage, associated with dense (10 4 cm"'), relatively 
cold (9.5 < Tk < 17 K ), and fairly quiescen t gas (typical line 
width s of 0.4 km s~\ iMuench etail 120071: IRathborne et alJ 
120081) . Non-thermal gas motions inside the cores are sub- 
sonic and mass independent. Therefore, thermal pressure ap- 
pears to be the dominant source of internal pressure. In addi- 
tion, these cores appear to be pres sure confined, but gravita- 
tionally unb ound (lLada et al.ll2008l) . 

Recently, lAlves et al.l (120081) performed an optical polari- 
metric survey toward the diffuse gas in the Pipe Nebula. They 
find a large scale magnetic field that appears to be mostly per- 
pendicular to the cloud main axis. The magnetic field ex- 
erts a pressure (~10 6 KcirT 3 ) that is likely responsible for 
driving the collapse of the gas and dust cloud along the field 
lines. The polarization properties significantly change along 
the Pipe Nebula. This fact allowed the authors to distinguish 
three regions in the cloud: B59, the stem, and the bowl (see 
Fig. Q]). B59 shows low polarization levels but high disper- 
sion of the polarization position angles. Moving through the 
stem toward the bowl, the polarization level increases and the 
dispersion decreases. These authors propose that these three 
regions might be in different evolutionary stages. B59 is the 
only magnetically supercritical region and the most evolved 
of the Pipe Nebula, the stem would be at an earlier evolution- 
ary stage, with material still collapsing, and finally, the bowl 
would be at the earliest stage, with cloud fragmentation just 
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Galactic Longitude 

Fig. 1. — Position of the observed co res plotted over the 2MASS extinction map of the Pipe Nebula (Lombardi et al, 2006). The dashed lines separate the three 
different magnetically defined regions I Alves et al. 2008). The lowest visual extinction ( Ay) corresponds to 0. 5 magnitudes. The highest A v is observed toward 
the bowl of the Pipe and the B59 region, where it reaches approximately 20 magnitudes (Lombardi et al. 2006). We selected cores located in all three regions of 
the Pipe (bowl, stem, and B59). 



TABLE 1 
Source list. 





a(J2000) 


<5(J2000) 


VLSR 




Source" 


h m s 


o I II 


(kms- 1 ) 


Region* 


Core 14 


17 12 34.0 


-27 21 16.2 


+3.6 


B59 


Core 40 


17 21 16.4 


-26 52 56.7 


+3.3 


Stem 


Core 48 


17 25 59.0 


-26 44 11.8 


+3.6 


Stem 


Core 109 


17 35 48.5 


-25 33 05.8 


+5.8 


Bowl 



(a) According to Lombardi et al. 1 2006) numbering. 

(b) According to lAlves et a l. (2008) diffuse gas polarimetric properties. 

started. 

Based on lAlves~eT al. (2008) results, we selected a sam- 
ple of cores distributed in the different regions of the Pipe 
Nebula. We started an extensive molecular survey of these 
cores using the IRAM 30-m telescope. The aim of this study 
is to probe their chemical evolutionary stage, which could 
be related with the dynam ical age according to chemical 
mode l ing of starless core s dTavlor et al.l 1 1998c iMorata et al.l 
120031: iTafalla et al.l 120061) . These models predict that some 
molecules, such as carbon-containing molecules, are formed 
very early in the chemical evolution, and are known as early- 



time molecules. These species are expected to be abun- 
dant in chemically young or low density cores, and most of 
the m seem to be affected earlier by depletion effects (see 
e.g.jTavlor et alJl998tlOhashi et alJl999ElBergin et al.ll2001b 
ITafalla et alJl2006h . Other species, such as nitrogen-bearing 
molecules and deuterated species, require a longer time to 
form. Thus, they are formed later in the chemical evolution 
and are known as late-time molecules. They are not expected 
to b e depleted until densities of 10 6 cm 3 are reached (see 
e.g. ICaselii et alj|2002t iFlower et allkOOd iBergin & Tafallal 
I2007t lAikawa et al.ll2008l)~ The qualitative comparison of the 
relative abundances of different types of molecules in each 
core can provide us with some clues about their possible evo- 
lutionary stage. From the observational point of view, there 
have been several authors that have studied the evolutionary 
stage of pre- and p rotostellar cores throu gh molecular sur- 
veys. For instance, iKontinen et al.l (120001) have used a large 
sample of molecules in a prestellar and a protostellar core. 
They find very different chemical compositions, specially in 
NaH + and long carbon-chain molecule abundances. The for- 
mer is typical of a pure gas-phase chemistry, while the lat- 
ter require an evolved chemistry to form. According to time- 
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dependent chemistry models they interpret the differences as 
differe nt stages of chemical evolution. Later, iTafalla et alj 
( 120041) made a chemical analysis of L1521E, which helped 
to determine the extreme yo uth of this prestellar core. From 
the theoretical point of view. lAikawa et all (12003) have simu- 
lated the evolution of a prestellar core and identified the dif- 
ferent molecular abundances at different evolut ionary stages 
to fina ll y compare the results with t he sample of Tafa lla et alj 
(120021) . iMorata et af] (I2003L 120051) have used the modeling 
results of iTaylor et al.l (119981) to compare with observations 
toward the L673 molecular cloud. 

Based on this, we observed a set of early- and late-time 
molecules (see Tabled toward the selected cores for a subse- 
quent comparison. In addition, we mapped the 1 .2 mm dust 
continuum emission of the cores to obtain a complete descrip- 
tion of the structure, chemistry and evolutionary stage of the 
four selected Pipe Nebula cores. 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. MAMBO-II observations 

We mapped the Cores 14, 40, 48, and 109 (according to 
the lLombardi et al1 l2006 numbering) at 1.2 mm with the 1 17- 
receiver MAMBO-II bolometer (array diameter of 240") of 
the 30-m IRAM telescope in Granada (Spain). The positions 
and velocity of the local standard of rest (vlsr) for each core 
are listed in Table Q] The observations were carried out in 
April and May 2009 and in January and March 2010, in the 
framework of a flexible observing pool. A total of 13 usable 
maps were selected for analysis: 3 for Cores 14, 40, and 109, 
and 4 for Core 48. The weather conditions were good, with 
zenith optical depths between 0. 1 and 0.3 for most of the time. 
The maps were taken at an elevation of <25° due to the decli- 
nation of the sources. 

The beam size of the telescope is ~11" at the effective fre- 
quency of 250 GHz. The sources were observed with the 
on-the-fly technique, with the secondary chopping between 
46" and 72" parallel to the scanning direction of the tele- 
scope. The telescope was constantly scanning at a speed of 
8" s for up to 65 s. This resulted in typical integration times 
for each map of ~1 hour. When possible, each source was 
mapped with different scanning directions (in equatorial co- 
ordinates) or rotating the secondary mirror of the telescope to 
avoid scanning artifacts in the final maps. We measured the 
zenith optical depth with a skydip and checked pointing and 
focus before and after each map. The average corrections for 
pointing and focus stayed below 3" and 0.2 mm, respectively. 
Flux density calibrators were observed every few hours. 

The data were reduced usin g MOPSIC with t he ite rative 
reduction strategy developed by Kauffmann et"ail(l2008l) . The 
main advantages of the new scheme are that (i) sources much 
larger than in the classical approach can be recovered, (ii) the 
signal-to-noise ratio (SNR) of the final map increases, and 
(Hi) they suffer from less artifacts. The figures were created 
using the GREG package, from the GILDASQ software. 

All the maps have been convolved with a 21. "5 Gaussian, 
larger than the telescope beam, in order to improve the SNR, 
and to smooth the appearance of the maps. The size of the 
Gaussian was chosen to be the one of the CN (1-0) molecular 
transition (see Tabled which provides good spatial resolution 
and large SNR for the four maps. 

1 MOPSIC and GILDAS data reduction packages are available at 
|http://www.iram.fr/'l RAMFR/GILDAS 



TABLE 2 

Molecular transitions observed toward the Pipe Nebula cores with the 
IRAM 30-m antenna. 







Frequency 


Beam 


Beam 


AV 




IVlOieCUie 


■ ■ 

1 1 iinsition 




y ) 


efficiency 


(km s" 1 ) 


type 


C,H, 


(2,,2-li,o) 


85.3389 


29.0 


0.78/0.81 


0.07 


E 


HCN 


(1-0) 


88.6318 


28.0 


0.78/- 


0.07 


E 


N 2 H + 


(1-0) 


93.1762 


26.5 


0.77/0.81 


0.06 


L 


C 34 S 


(2-1) 


96.4130 


26.0 


-/0.81 


0.06 


E 


cs 


(2-1) 


97.9809 


25.5 


0.76/0.81 


0.06 


E 


CN 


(1-0) 


113.4909 


21.5 


0.75/0.81 


0.05 


E 


N 2 D + 


(2-1) 


154.2170 


15.0 


0.77/0.74 


0.04 


L 


DCO+ 


(3-2) 


216.1126 


10.5 


0.57/0.63 


0.03 


L 


CN 


(2-1) 


226.8747 


10.0 


0.53/0.63 


0.03 


E 


N 2 D + 


(3-2) 


231.3216 


10.0 


0.67/0.63 


0.03 


L 


H"CO + 


(3-2) 


260.2554 


9.0 


0.53/0.63 


0.02 


L 


(a) ABCD and EMIR receiver, respectively 

(b) Spectral resolution. 



(c) E = Early-time. L — Late-time. See i^and 34.3l for details. 

2.2. Line observations 

We made several pointed observations within the regions 
of the Cores 14, 40, 48, and 109 with the heterodyne re- 
ceivers of the 30-m IRAM telescope (ABCD and EMIR re- 
ceivers). The observations were carried out in three epochs. 
The first epoch was August and September 2008. We used 
the capability of the telescope to perform simultaneous ob- 
servations at different frequencies to observe the emission 
of the C 3 H 2 (2i 2-1 1 0), HCN (1-0), N 2 H+ (1-0), CS (2-1), 
CN (1-0), N 2 D+ (2-1), DCO + (3-2), CN (2-1), N 2 D+ (3-2) 
and H 13 CO + (3-2) molecular transitions arranged in 3 dif- 
ferent frequency setups covering the 3, 2, 1.3 and 1.1 mm 
bands. To do this, we combined the A100/B100/A230/B230 
and A100/D150/A230/D270 SIS heterodyne receivers. The 
observational strategy was first to observe several positions 
with a 20 " spacing cent e red on the C ls O pointing center re- 
ported by iMuench et alj (120071) (depicted by star symbols in 
Fig.|2]i, which is v ery close to the visual extinction peak po- 
sition of each core dLombardi et alj|2~006t) . The visual extinc- 
tion peak is assumed to be the de nsest region of the cor e, and 
it was defined as core center by IMuench et alj d2007l) . The 
second and third epochs were August 2009 and June 2010, 
respectively, both using the new EMIR E0/E1/E2 receivers. 
We observed deeper toward the position of the grid of first 
epoch closer to the dust continuum peak (see circle symbols 
in Fig. [2j». We observed also the C 34 S (2-1) molecular transi- 
tion. Table [2] shows the transitions and frequencies observed. 
We used the VESPA autocorrelator as the spectral backend, 
selecting a channel resolution of 20 kHz, which provided a 
total bandwidth of 40 MHz. The corresponding velocity res- 
olutions, main-beam efficiencies and half-power beam widths 
at all the observed frequencies are also listed in Table [2] We 
used the frequency-switching mode with a frequency throw 
between 3.83 and 22.98 MHz, depending on the transition. 
System temperatures in nights considered "good" were be- 
tween 200 to 275 K at 3 mm, and between 440 and 960 K at 
1 mm (T sys reached 450 K and 3200 K in bad nights, respec- 
tively). Pointing was checked every 2 hours. 

We reduced the data using the CLASS package of the 
GILDAS 1 software. We obtained the line parameters either 
from a Gaussian fit or from calculating their statistical mo- 
ments when the profile was not Gaussian. 

3. RESULTS AND ANALYSIS 
3.1. Dust continuum emission 
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Fig. 2. — IRAM 30-m MAMBO-II maps of the dust continuum emission at 1 .2 mm toward four cores of the Pipe Nebula. The greyscale levels for all the maps 
are 3 to 18 times 5.75 mJy beam -1 . The contour levels are 3 to 1 1 times a in steps of 1-cr, for Cores 14, 40 and 48, and 3 to 21 cr in steps of 3-tr for Core 109. 
1-cr is 4.5, 5.0, 3.5, and 4.5 mJy beam -1 for Core 14, 40, 48, and 109, respectively. The red thick contour mark the half maximum emission level of the source 
(see Table[3) . Black or white stars indicate the C ls O pointing center reported by Muench et al. ( 2007), which is very close to the visual extinction peak position 
of each core iLombardi et al. 2006). Black or white filled circles indicate the position where line observations have been performed , close to the dust continuum 
emission maximum which falls into the beam area. The blue vectors depict the magnetic field vector found by Franco et al. 1 2010). Note that for Core 40 there 
are no optical polarimetry measurements on the east side due to the high visual extinction. In the bottom left corner of the bottom right panel the convolved beam 
and the spatial scale for the maps is shown. 



TABLE 3 

1.2 MM CONTINUUM EMISSION PARAMETERS. 





a(J2000) " 


<S(J2000) 


Tdust 


RMS 


Sy 


^Pcak 


Diameter 


N» 2 " 


nn 2 b 


Mass b 


Source 


h m s 




(K) 


(mJy beam -1 ) 


(Jy) 


(mJy beam -1 ) 


(PC) 


(10 21 cm -2 ) 


(10 4 cm -3 ) 


(M Q ) 


Core 14 (filament) 


17 12 31.5 


-27 21 41.0 


12.0 C 


4.5 


2.56 


51.6 


0.106 


12.21 


5.59 


2.87 


Core 14 (core) 










1.24 




0.071 


13.27 


9.09 


1.40 


Core 40 


17 21 14.7 


-26 52 47.8 


10.3'' 


5.0 


1.73 


42.0 


0.104 


11.05 


5.16 


2.51 


Core 48 


17 25 57.3 


-26 44 22.3 


10.0 rf 


3.5 


1.44 


27.9 


0.127 


6.14* 


2.35 rf 


2.09 d 


Core 109 


17 35 47.7 


-25 32 52.9 


9.5 C 


4.5 


2.76 


105.3 


0.063 


47.60 


36.57 


4.00 



(a) Pointing position of the chemical observations which lies inside the same beam area of the dust continuum emission peak. 

(b) Assuming K250 g//j=0.0066 cm 2 g -1 as a medium value between dust grains with thin and thick ice mantles I Ossenkopf & Hennina 1994). See Appendix IaI 
for details on the calculation. 

(c) Adopted to be equal to the kinetic temperature derived for NH3 (Rathborne et al. 2 0081) . 

(d) No kinetic temperature estimate, therefore we assumed 10 K based on the temperatures of the other cores IRathborne et al.120081) . 
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TABLE 4 

H2 COLUMN DENSITIES", A^Ht , OF THE PlPE NEBULA CORES IN CM~ 



Source 


10." 5 


15."0 


21. "5 


27."0 


Core 14 


1.75 x 10 22 


1.38 x 10 22 


1.21 x 10 22 


1.11 x 10 22 


Core 40 


1.32 x 10 22 


1.28 x 10 22 


1.12 x 10 22 


1.07 x 10 22 


Core 48 


1.09 x 10 22 


8.88 x 10 21 


7.38 x 10 21 


6.99 x 10 21 


Core 109 


4.19 x 10 22 


3.73 x 10 22 


3.23 x 10 22 


3.08 x 10 22 



(a) Average column densities are calculated within one beam area toward 
the dust continuum emission peak. The values of K750 GHz ar >d ^dust ale me 
same as for Table [3] These values are combined with the molecular column 
densities to find the molecular abundances in the same beam area. The corre- 
spondence is: 10." 5 with DCO+, 15."0 with N 2 D + , 21." 5 with CN (1-0) and, 
finally, 27."0 with C 3 H 2 , HCN, N 2 H + , CS and C 34 S. 

In Fig. |2] we present the MAMBO-II maps of the dust con- 
tinuum emission at 1 .2 mm toward the four selected cores of 
the Pipe Nebula, convolved to a 21. "5 beam. Table[3]gives the 
peak position of the 1.2 mm emis sion after convolution with 
a Gaussian, the dust temperature dRathborne et alj|2008l) . the 
RMS noise of the emission, the flux density and the value 
of the emission peak. Additionally, we also give the derived 
FWHM equivalent diameter, which is the diameter of the cir- 
cular area equal to the area within the FWHM level, depicted 
by a red contour in Fig. [2] Table [3] also lists the H2 column 
and volume density, as well as the mass for each core. These 
parameters are derived from the emission within the 3-<x level 
and discussed in ^4] 

The flux density of the cores ranges between ~1.24 and 
~2.76 Jy. Note, however, that the extinction maps show that 
the studied cores are surrounded by a diffuse medium (see 
Fig. [Tl and [Lombardi et al.l 12006b . The on-the-fly reduction 
algorithms assume that the map limits have a zero emission 
level. Due to the presence of the diffuse material, this could 
not be true for the observed cores, and, therefore, the mea- 
sured flux density of the maps might be lower than the actual 
value. We derived average H2 column densities (A^h,, see Ap- 
pendix [A} toward the dust continuum emission peak for the 
different resolutions (listed in Table |4]i of the detected molec- 
ular transitions (see Table |5). We derived their abundances 
with respect to H2. The results are discussed in §|4] 

The maps of Fig. [2] show th e different morpho logy of the 
cores. Following the results of lAlves et al.l (120081) . it is inter- 
esting to compare the shape of the cores with their location 
along the Pipe Nebula. Core 14, located in B59, belongs to a 
clumpy and filamentary structure of ~500" (~0.35 pc) elon- 
gated along the NE-SW direction. T his is in perfect agree- 
ment with previous extinc tion maps (lLombardi et al.l l2006t 
iRoman-Zuniga et al.ll2009h . On the other hand, Core 109, lo- 
cated in the bowl, shows a compact and circular morphology 
with a FWHM of -90" (-0.063 pc). Cores 40 and 48, both lo- 
cated in the stem, have elliptical morphologies with extended 
diffuse emission. 

3.2. Molecular survey of high density tracers 

iMuench et al.l d2007l) reported C ls O pointed observations 
toward the Pipe Nebula cores measured with a resolution of 
56". As seen in their Fig. 1, the position of the C ls O is 
very close to the visua l extinction peak position of each core 
( lLombardi et al.1120061) . Our higher resolution maps show a 
peak position offset for all the cores. As listed in Tables Q] 
and [3] and as shown in Fig. [2] the dust continuum peak does 
not coincide exactly with the Ay peak (stars in Fig. |2j. How- 
ever, the difference is compatible with the angular resolution 
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Velocity (km/s) 

Fig. 3. — IRAM 30-m line spectra of the molecular transitions with no 
hyperfine components toward the four selected cores of the Pipe Nebula (see 
Table [T}. The name of the core is indicated above the top panel of each 
column. Rows correspond to a single molecular transition specified on the 
second column. The velocity range is 5 km s , and it is centered on the 
vlsr of each core is marked with a vertical dotted line. Vertical axis show 
the 7mb of the emission, and the zero level is marked by a horizontal dotted 
line. Some of the Core 109 spectra, with the highest measured 7/mb, have 
been divided by 2 to fit to the common scale. 

of the extinction maps. We decided to present only molec- 
ular line data of the observed positions closer to the dust 
continuum emission peak (circles in Fig. |2]l , defined as the 
core center and supposed to exhibit brighter emission from 
molecular transitions. The typical core size is ~90" or larger 
(see Table [3). The beam size of the detected lines, except 
for N 2 D + (2-1) and DCO + (3-2), range from 21."5 to 29."0, 
while the initial grid of pointed position had a separation of 
20."0, thus the emission peak stays within the beam area for 
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Fig. 4. — IRAM 30-m line spectra of the molecular transitions with hyperfine components toward the four selected cores of the Pipe Nebula (see TablefT). The 
name of the core is indicated above the top panel of each column. Rows correspond to a single molecular transition specified on the second column. The velocity 
range is 16.5, 20 and 12 km s _1 for HCN (1-0), N 2 H+ (1-0) and CN (1-0) respectively. The vlsr of each core is marked with a vertical dotted line. Vertical axis 
show the Tmb of the emission, and the zero level is marked by a horizontal dotted line. Core 109 spectra, with the highest measured Tub, have been divided by 
2 to fit to the common scale. 



TABLE 5 

Summary of detections and upper levels in K toward the Pipe Nebula 
cores" . 



Molecular 



Core 



transitions 


14 


40 


48 


109 


C 3 H 2 (2i,2-li, ) 


V 


V 


<0.07 


V 


HCN (1-0) 


<0.21 


V 


V 


V 


N 2 H + (1-0) 


V 


V 


<0.07 


V 


C 34 S (2-1) 


V 


V 


V 


V 


CS (2-1) 


V 


V 


V 


V 


CN (1-0) 


V 


V 


<0.17 


V 


N 2 D + (2-1) 


<0.12 


V 


<0.08 


V 


DCO+ (3-2) 


<1.71 


<0.61 


<0.76 


V 


CN (2-1) 


<0.97 


<1.70 


<0.76 


<0.90 


N 2 D + (3-2) 


<1.01 


<0.93 


<1.94 


<0.91 


H 13 CO + (3-2) 


<1.52 


<1.40 


<2.38 


<1.34 



molecular transitions at 3 mm (C3H2 (2i j2 -l 1,0), HCN (1-0), 
N 2 H + (1-0), C 34 S (2-1), CS (2-1), and CN (1-0)), although 
their intensities are lower than for Core 109. Core 14 shows 
emission in all the 3 mm transitions except in HCN (1-0). Fi- 
nally, Core 48 only shows emission in CS (2-1), C 34 S (2-1), 
and HCN (1-0). 

In addition to the line parameters, we derived the molecular 
column densities for all the detected species (see AppendixlBl 
for details) which are listed in Table [7] For the transitions 
with detected hyperfine components (HCN, N 2 H + , and CN) 
we derived the opacity using the hyperfine components fitting 
method of the CLASS package. For the CS and C 34 S molec- 
ular transitions we derived numerically the opacity using 



r MB (C M S) 1 - exp(-r) 



(a) The transitions marked with y have been detected toward the correspond- 
ing core. Otherwise, the 3<x upper limit is shown. 

these molecular transitions. Therefore, the molecular line 
properties that we obtain are representative of the chemistry 
of the core center. 

Table|5]summarizes the detections or the 3cr upper limits of 
the non detections toward each core. Table[6]give the parame- 
ters of the detected lines. In Figs.[3]and|4]we show the spectra 
of the different molecular transitions observed toward the dust 
continuum emission peak of each core. Core 109 shows the 
stronger emission in all the detected transitions in our sam- 
ple. This is the core with the most compact and circular mor- 
phology (see Fig. |2]). Core 40 shows also emission in the six 



r MB (CS) 1 - exp(-r r) ' 



(1) 



where r is the CS to C 34 S abunda nce ratio, assumed to be 
equal to the terrestrial value (22.5. iKim & Kodl2003h . We 
found a high opacity toward Cores 14 and 48 for CS (2-1), 
10.8 and 6.0 respectively, whose spectra show self-absorption 
(see Fig. [3]). For Cores 40 and 109 we found lower opaci- 
ties, r=3.1 and 4.2 for CS (2-1) respectively. We assumed 
optically thin emission in C 3 H 2 , DCO + , and N 2 D + (2-1), the 
latter with only the main hyperfine component detected. This 
conservative assumption could not be true, so the column den- 
sities should be taken as lower limits. We also derived the 
molecular abundances with respect to H 2 (see Table [8}, tak- 
ing into account the resolution for each molecular transition 
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TABLE 6 

Line parameters . 



Molecular 




Tub b 


A xt c 


J'TMBdv* 


vlsr 


Av LS r 






transition 


Source 


(K) 


(K) 


(K kms" 1 ) 


(km s" 1 ) 


(kms -1 ) 




Profile'' 


C3H9 (2j 2 — 1 1 0) 


Core 14 


0.37(6) 




0.086(11) 


3.502(14) 


0.22(3) 




G 




Core 40 


1.19(5) 




0.347(9) 


3.420(4) 


0.273(9) 




G 




Core 109 


2.74(6) 




0.799(13) 


5.8340(20) 


0.274(5) 




G 


CS (2-1) 


Core 14 


0.69(10) 




0.41(3) 


3.439(21) 


0.45(4) 


10.8(1.1) 


SA 




Core 40 


1.94(7) 


_ 


0.560(17) 


3.369(4) 


0.415(14) 


3.1(3) 


NS 




Core 48 


0.79(7) 




0.402(18) 


3.684(11) 


0.477(22) 


6.0(6) 


SA 




Core 109 


1.93(8) 


_ 


0.743(17) 


5.836(4) 


0.361(9) 


4.2(4) 


G 


C 34 S (2-1) 


Core 14 


0.267(25) 




0.068(5) 


3.545(8) 


0.241(20) 


0.5(1) 


G 




Core 40 


0.268(16) 




0.069(3) 


3.381(5) 


0.241(13) 


0.14(1) 


G 




Core 48 


0.187(23) 


_ 


0.041(4) 


3.729(11) 


0.20(3) 


0.26(3) 


G 




Core 109 


0.34(3) 




0.083(5) 


5.825(7) 


0.233(17) 


0.19(2) 


G 


N 2 D+ (2-1) f 


Core 40 


0.084(20) 




0.019(3) 


3.280(15) 


0.21(3) 




G 




Core 109 


0.31(4) 




0.109(7) 


5.673(11) 


0.331(22) 




G 


DCO+ (3-2) 


Core 109 


0.70(11) 




0.151(18) 


5.828(13) 


0.202(21) 




G 


HCN(l-O) 


Core 40 




1.55(11) 




3.410(16) 


0.334(22) 


6.0(5) 


NS 




Core 48 


_ 


0.33(10) 




3.54(5) 


0.90(11) 


2.4(1.2) 


G 




Core 109 (1) 




2.53(3) 




5.93(7) 


0.16(22) 


0.25(10) 


NS 




Core 109 (2) 




6.10(3) 




5.72(7) 


0.25(22) 


10.20(10) 


NS 


N 2 H + (1-0) 


Core 14 




0.0341(16) 




11.500(5) 


0.206(10) 


0.10(9) 


G 




Core 40 




0.219(12) 




11.4000(19) 


0.249(5) 


0.171(25) 


G 




Core 109 




0.904(14) 




13.8000(5) 


0.2150(11) 


0.467(11) 


G 


CN(l-O) 


Core 14 




0.051(9) 




3.64(8) 


0.81(15) 


0.1(7) 


G 




Core 40 




0.65(22) 




3.430(21) 


0.36(5) 


3.9(1.3) 


G 




Core 109 (1) 




1.41(22) 




5.930(5) 


0.162(11) 


1.13(23) 


G 




Core 109 (2) 




2.3(1.3) 




5.670(7) 


0.101(16) 


4.(3) 


G 



(a) Line parameters of the detected lines. The former five molecular transitions have no hyperfine components. The parameters for the transitions labeled as 
G (see last column) have been derived from a Gaussian fit while line parameters of NS and SA profiles have been derived from the intensity peak (Tub), and 
zero (integrated intensity), first (line velocity) and second (line width) order moments of the emission. The latter three molecular transitions have hyperfine 
components. The parameters have been derived using the hyperfine component fitting method of the CLASS package. The value in parenthesis shows the 
uncertainty of the last digit/s. If the two first significative digits of the error are smaller than 25, two digits are given to better constrain it. 

(b) Only for molecular transitions with no hyperfine components. 

(c) Only for molecular transitions with hyperfine components. 

(d) Derived from a CLASS hyperfine fit for molecular transitions with hyperfine components. Derived numerically for CS and C 34 S using Eq.fT] A value of 0.3 
is assumed when no measurement is available. 

(e) G: Gaussian profile. NS: Non-symmetric profile. SA: Self-absorption profile. 

(f) Only the main component is detected. 



(see Table©. 

4. DISCUSSION 

We observed four selected cores located in the different re- 
gions of the Pipe Nebula (bowl, stem, and B59), in different 
molecular tracers and in dust continuum emission, to study 
and compare their physical and chemical properties. The 
cores were selected based on th e results of the op tical po- 
larimetric survey carried out by lAlves et afl (120081) . In the 
following subsections, we discuss and compare the proper- 
ties of each individual core, as well as an overall analysis of 
such properties, and we try to relate our results with previous 
works. In particular, in the next subsection we compare the 
dust continuum emissi on with the visual extinction maps of 
IRathborne et alJ (120081) and the trend found for the diffuse gas 
:al.l (120081) . 



by Al 



ves et ; 



4.1. Comparison of visual extinction and 1.2 mm continuum 
emission maps 

The beam size of our observations is 1 1", c onvolved to a 
Gaussian of 21. "5 in th e maps shown (see §13. lb . while that of 
IRathborne et alJ (12008b is ~60". Hence, our maps suffer from 
less beam dilution and we can resolve smaller structures. The 
sensitivity limit o f the Ay observations is fixed at 1.2 mag 
dLada et alj|2008l) . which corresponds to a column density of 



~lxl0 cm (IWagenblast & Hartquistiri989h . A conserva- 
tive estimation of the sensitivity limit of our maps, in the same 
conditions, can be derived using the 3-cr emission level of the 
noisiest continuum map convolved to a 60" Gaussian. The 
resulting beam averaged column density, for a Tk of 10 K, is 
~4xl0 cirT 2 . Therefore, as seen from the minimum column 
densities in th e same conditions, our dataset have slightly bet- 
ter sensitivity. lLada et aTl ( 12008b defines the equivalent radius 
of the core using the region with emission brighter than 3-cr, 
while we use the region with emission brighter than half of 
the peak value. This difference prevents a direct comparison 
of the radii and densities. The core masses, however, depend 
only on the integrated flux density and can be compared. Our 
masses are on average ~3.4 times smaller, ranging from ~0.9 
for Core 109 to ~7 for Core 14. 

We esti mated the difference b etween the 2MASS extinc- 
tion maps (lLombardi et al.ll2~006l) and the 1.2 mm dust maps. 
To do this we first transformed the original near-IR extinc- 
tion maps to visua l extinction maps using Ay = Ajj/0.118 
jDutraetal.ll2002h . Then, we convolved the 1.2 mm dust 
maps with a Gaussian of 60" to have the same resolution. 
We transformed the 1.2 mm dust maps to column density 
maps (see Appendix lAl for details). We assumed a uniform 
typical temperature of 10 K for all the cores. To estimate 
the uncertainty caused by this assumption, we also made 
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-2.65 -2.7 -2.75 -2.65 -2.7 -2.75 



Fig. 5. — Left panels: C olor image of the visua l extinction map derived from 
the near-IR observations (Lombardi et al. 2006) superposed with the contour 
map of the extinction map derived from our 1.2 mm dust continuum maps in 
galactic coordinates. Contours are from 2.5 to 30 visual magnitudes by steps 
of 2.5. Right panels: Color image of the difference between the visual extinc- 
tion map derived from the near-IR and the 1.2 mm observations within the 
region of the A v converted dust continuum maps with A v >2.5 magnitudes. 
Contours are the same as in the left panels. Core number is indicated in the 
lower right corner of the panels. The color scale (in visual magnitudes) are 
shown in the right side of the panels. 

the calculations for temperatures of 8 and 12 K, which re- 
sulted in an average maximum variation over the whole map 
of ~2.4 extinction magnitudes. We also assumed for all the 
cores #250 gh<-=0.0066 cm 2 g~' as average value between dust 
grains with thin and thick ice mantles for a volume density 
of ~10 5 cm 3 (lOssenko pf & Hennind fl994) with an uncer- 
tainty of about a factor of 2. As a final step, we used the 
relatio nship Ay = 1.258 x 1(T 2i ./Vh 2 (IWagenblast & Hartquisfl 
119891) to transform column density to visual extinction. The 
resulting maps of the difference between the extinction maps 
derived from near-IR and mm data are shown in the right- 
hand side panels of Fig. [5] For Core 40 we found an ex- 
cess of extincti on t hat could be due to the filtered diffuse 
emission (see 33. 11 1. However, in such a case, one would 
expect this excess to be present over the whole map. For 
Cores 14 and 48 we found a good agreement between both 
tracers. On the other hand, at denser regions such as the cen- 
ter of Core 109 («h 2 >4x10 5 cm -3 ), the Ay derived from the 
1.2 mm dust is significantly larger, >10 mag, than that derived 
from the near-IR. This is the core with the highest column 
density (~4.8xl0 22 cm -2 ), therefore this suggests that near- 



IR extinction maps constructed from 2MASS catalogs do not 
have enough sensitivity or sampling scale to resolve the cen- 
ters of very dense cores. In such dense regions, the num- 
ber of 2MASS catalog background stars is not high enough 
to provide neither a large number of sources per pixel, nor 
a large number of high extinction measurements, thus the 
high extinction regions might be poorly resolved and under- 
estimated. These biases may explain, co mbined with the 
larger radius, the lower densities reported by Rat hborne et al.l 
(120081) . Extensive observations toward the Perseus cloud in 
visua l extinction and i n radio continuum provide similar re- 
sults dKirk et al.ll2006l) . Extinction maps with higher resolu- 
tion, made with deeper observations, are able to resolve better 
the high extinction levels. For example. iKandori et al.l (120051) 
observed Core 109 (named FeSt 1-457) in Ay deeper with a 
resolution of ~30", and found a morphology in perfect agree- 
ment with our continuum observations. Their Ay intensity 
peak at the core center of Ay~41.0 (the largest in their sam- 
ple) is very close to our derivation, Ay~39.2, for a 30" beam. 

Summarizing, our dust continuum maps seem to be better 
at tracing the high extinction regions of the prestellar cores, 
at least at this spacial resolution. These results suggest that 
the dust continuum emission would trace the dense and cold 
cores better than the 2MASS derived visual extinction. On 
the other hand, the visual extinction would be more sensitive 
to the cloud diffuse extended emission. 

4.2. Discussion on the individual cores 
4.2 A. Core 14 

Core 14, located in B59, is a compact and dense core but the 
less massive in our sample. It is the only core that belongs to 
a clumpy and filamentary structure, which is elongated along 
the NE-SW direction with an extent of ~500" (~0.35 pc, see 
Fig. |2]), with a morphology quite similar to that shown in Ay 
maps lRathborne et al.l2008blRoman-Zuniga et alJ2009l) . The 
location of Core 14 inside an elongated and clumpy filament, 
suggests that probably it is still undergoing fragmentation, 
which could lead to the formation of smaller cores. In fact, 
it is resolved in several small clumps which have sizes com- 
parable to the sizes of the other cores in the Pipe, with radii of 
about ~0.04 pc. 

Core 14 shows emission in all the early-time molecules at 
3 mm. CS (2-1) and C 34 S (2-1) are clearly detected (see 
Fig- Eb-. an d the abundances are the largest of the sample (see 
Table©. On the other hand, C 3 H 2 (2 u -li,o) and N 2 H + (1- 
0) show weak emission and, consequently, low abundances. 
iRathborne et alj (12008b detect weak NH3 emission, in good 
agreement with ourN 2 H + measurements. Only the main com- 
ponent of the CN (1-0) transitions is clearly detected. These 
are signatures of an object very young chemically. 

4.2.2. Core 40 

Core 40, located in the stem, is another core with irregular 
morphology. This core shows emission in all the transitions at 
3 mm, of both early- and late-time molecules, and in the late- 
time N 2 H + transition at 3 mm. The molecular emission of 
Core 40 is strong, and only the emission of Core 109 is more 
intense, except for CS, which shows the same Tmb for both 
cores. These cores are the only ones that show strong late- 
time molecule emission. Core 40 presents the highest CN and 
N 2 H + abundances (see Table©. Regarding N 2 H + , the intense 
emission with all the hyperfine components d etected is in per- 
fect agreement with previous results of NH3 ( Rat hborne et alj 
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[2008b . The HCN emission for Core 40 is quite anomalous, be- 
cause the main hyperfine component is weaker than the satel- 
lite compone nts (see Fig.B}. This suggests that the em ission 
is not in LTE. iGonzalez- Alfonso & Cernichard (119931) inves- 
tigated with Monte-Carlo techniques the variation in HCN (1- 

0) profiles. According to their work, an infalling cloud with 
a dense c entral core (see Fig. [Zjl surrounded by a large diffuse 
envelope (lLombardi et al.|[2006b may produce an HCN (1-0) 
spectrum as the observed toward Core 40. 

4.2.3. Core 48 

Core 48, located in the stem, has a quite elongated mor- 
phology. It is embedded in an environmen t with high polar- 
ization angle dispersion dAlves et al.M2.Q08l) . which is the ex- 
ception of this polarimetrically defined region. It is very dif- 
fuse, this is the largest and the less dense core in the sample. 
It shows emission only in three early-time molecules: CS (2- 

1) , C 34 S (2-1) and, marginally, HCN (1-0). The abundances 
of CS and C 34 S are among the largest in the sample, slightly 
lower than those for Core 14. The N2H + molecule was un- 
detected, in agreement with previous measurements of NH3 
dRathborne et al.ll2008b . 

4.2.4. Core 109 

Core 109, located in the bowl, is the most circular and com- 
pact core in the sample. The dust continuum emission of this 
core is similar to that of the other cores. However, it is the 
densest one in our s ample and the most massive (~4 M ). 
iKandori eT al. (2005) find, through a Bonnor-E bert profile 
fit, th at Core 109 is gravitationally unstable. lAguti et all 
(12007b find, through observations of molecular transitions, 
that thi s core (designated als o as FeSt 1-457) is gravitationally 
bound. IKandori et al.l (120051) suggest other models, apart from 
a Bonnor-Ebert sphere, including extra supporting mecha- 
nisms that might fit the density profile. lAguti et all (12007b 
propose that Core 109 is pulsating, based on expansion mo- 
tions of the outer layers. However, their Jeans mass measure- 
ment is compatible with the mass of the core and they propose 
a quasi-stable state near hydrodynamic equili brium . This core 
is embedded in a magnetized medium (see § 14.41) . thus, mag- 
netic support could be a plausible source of external support. 

This core shows emission of all the detected early- (C3H2, 
HCN, CS, C 34 S, and CN) and late-time molecules (N 2 H + , 
N 2 D + and DCO + ). The molecular emission of this core 
is always the strongest. Core 109 shows a very strong 
N 2 H + emission, i n agre ement with the NH3 measurements by 
IRathborne etail (12008b . As seen in Table [8] Core 109 has 
similar abundances for early-time molecules to those of the 
other cores. Interestingly, the CS and C 34 S abundances are 
the lowest in our sample, whi ch suggests CS dep letion toward 
the center (detected on C 18 Q. lAguti et al.ll2007l) . 

4.3. Qualitative chemistry analysis 

Table |4] shows a variation of about a factor of ~4 around 
10 22 cirT 2 of the average H 2 column densities derived for 
each of the cores with a 27" beam, the one used to calcu- 
late the abundances for the molecular transitions at 3 mm. 
This represents, using the relati onship Ay = 6.289 x 10~ 22 Nh 
(IWagenblast & Hartquistl fl989). average values of Ay~4.4 to 
~19.4. The first case would represent a shallow core, more 
affected by the external radiation field, which tends to have 
a younger chemistry. The other extreme probably indicates 
a denser and more shielded core, where one would expect to 



find more complex and evolved molecules. However, note 
that t his also depends on the time-scale n eeded to form the 
core dTafalla et alJl200l iCrapsi et alJl2005h . 

We find that CS (see Table [H), an early-time molecule, is 
detected in all the cores with abundances with respect to H 2 
of a fe w times 10" , sim ilar to the ones found in other dense 
cores (llrvine et"al]|1987b or the ones obtained in gas- phase 
chemical models dTavlor et al.lll998H"Garrod et alj|2004b . It is 
worth to mention that Cores 14 and 48 show high CS abun- 
dance, one order of magnitude higher than Cores 109 and 40. 
A similar result is found for the C 34 S abundances. The de- 
rived abundances for the early-time molecule HCN toward the 
cores in our sample is very uniform, and seems to be inde- 
pendent of their physical properties. The early-time molecule 
CN, a molecule that is also detected commonly in dense cores, 
has also a significantly lower abundance (a factor >4) toward 
Core 48 than toward the rest of the sample. Where detected, 
the CN abundance varies only within a factor of five. On the 
other hand, another early-time molecule such as C3H 2 , shows 
differences in abundances of at least a factor of five among 
Cores 14 and 48 with respect to Cores 40 and 109. Late-time 
molecules, such as N 2 H + or deuterated molecules, are not 
broadly detected in our sample: N 2 H + is detected except in 
Core 48. On the contrary, N 2 D + is only detected on Cores 40 
and 109, and DCO + only in Core 109. We found a higher 
abundance of N 2 H + t oward Core 40 than tow ard Core 109 by 
a factor of ~2, while IRathborne et al.l (120081) found an abun- 
dance of NH3 toward Core 109 higher than that of Core 40 by 
a factor of ~3.4. However, both Cores 40 and 109 show higher 
abundances in N 2 H + than Cores 14 and 48. Despite Cores 14 
and 40 having a similar average column density, the former 
shows 5 times less abundance of N 2 H + than the latter. More- 
over, Core 14 does not show emission in any other late-time 
molecule while Core 40 is detected in N 2 D + showing an abun- 
dance only a factor of 4 lower than that of Core 109. Briefly, 
the higher abundances in Cores 109 and 40 with respect to 
Core 14 and in particular to Core 48 (except for CS) is an in- 
dication that Cores 109 and 40 are more chemically evolved 
than Cores 14 and 48. However, the molecular abundances of 
these two late-time species are roughly an order of magnitude 
lower than the proto typical starless cores L1517B and L1498 
dTaf alia etal .112006b . which suggests that Cores 109 and 40 
m ay be in an earlie r evolu tionary stage than cores in Taurus. 

IRathborne et al.l d2008b observed the emission of the 
NH 3 (1,1), NH 3 (2,2), CCS (2i-l ), and HC 5 N (9-8) tran- 
sitions towards 46 cores of the Pipe Nebula. Cores 14, 40, 
48, and 109 were included in their observations. None of 
the lines were detected in Core 48, which is shown to be 
again the more chemically poor core of our sample. HC5N 
was not detected in core 14, which also has the weakest CCS 
and NH3 lines. The four transitions were detected in cores 40 
and 109, but with some differences. The NH3 lines are much 
more intense in core 109, a factor of ~4 for the (1,1) tran- 
sition and ~9 for the (2,2) line, while the CCS line is more 
intense in core 40, less than a factor of ~2, and the HC5N 
lines are very similar in both cores, inside the RMS. All these 
results are consistent with our observations: Core 48, which 
did not show emission of late-time molecules, is very poor 
chemically and shows a very young chemistry. Core 14, has 
some very weak emission of late-time molecules (NH3) but 
only weak emission of CCS, an early-time molecule. Core 
40 is more evolved chemically and shows stronger emission 
of early-time molecules than of late-time molecules. Finally, 
Core 109 is the one showing more diversity of molecules 
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TABLE 7 

Molecular column densities of the chemical species observed toward the Pipe Nebula cores fn cm~ 



Source 


C 3 H 2 " 


CS 


c 34 s 


CN 


HCN 


N 2 H + 


N 2 D +fl 


DCO+ " 


Core 14 


3.84 X 10" 


3.07 x 10 13 


6.19 x 10" 


1.16 x 10 12 


< 5.42 x 10 10 


9.70 x 10 10 


< 8.09 x 10 09 


<5.13x 10" 


Core 40 


1.65 x 10 12 


7.11 x 10 12 


2.94 x 10" 


4.70 x 10 12 


2.57 x 10 12 


4.89 x 10" 


2.17 x 10 09 


<4.93x 10 10 


Core 48 


<7.93x 10 10 


1.53 x 10 13 


2.94 x 10" 


< 1.64 x 10" 


2.59 x 10 12 


< 3.79 X 10'° 


< 6.88 x 10 09 


< 7.82 X 10 10 


Core 109 


6.36 x 10 12 


1.25 x 10 13 


3.68 x 10" 


2.76 x 10 12 


9.80 x 10 12 


6.79 x 10" 


2.72 x 10 10 


1.41 X 10" 



(a) Transition with no opacity mesurements available, thus optically thin emission is assumed to obtain lower limits of the column densities. 

TABLE 8 

Abundances of the chemical species with respect to H 2 observed toward the Pipe Nebula cores. 



Source 


C 3 H 2 b 




CS 




C 34 S 




CN 




HCN 




N 2 H + 




N 2 D + " 


DCO+ a 




Core 14 


3.45 x 10" 


ii 


2.77 X 10 


-09 


5.58 x 10- 


ii 


9.58 x 10" 


ii 


< 4.88 x 10" 


-12 


8.73 X 10" 


-12 


< 5.84 x 10~ 13 


< 2.93 X 10" 


ii 


Core 40 


1.55 x 10- 


■10 


6.64 x 10' 


-10 


2.75 x 10- 


-ii 


4.19 x 10" 


-10 


2.41 x 10" 


-10 


4.58 x 10" 


-11 


1.69 x 10~ 13 


< 3.73 X 10" 


-12 


Core 48 


< 1.13 x 10" 


-ii 


2.19 x 10" 


-OS 


4.21 x 10- 


-ii 


< 2.22 x 10" 


-11 


3.71 x 10" 


-10 


< 5.43 x 10" 


■12 


<7.75x 10~ 13 


<7.17x 10" 


-12 


Core 109 


2.06 x lQ- 


■10 


4.06 x 10" 


-10 


1.19x lQ- 


•11 


8.54 x lQ- 


-11 


3.18 x 10" 


■10 


2.20 x 10" 


-11 


7.30 x 10~ 13 


3.37 x lQ- 


■12 



(b) Transition with no opacity mesurements available, thus optically thin emission is assumed to estimate a lower limit of the column densities and, consequently, 
of the abundances. 

(c) Due to the lack of C 34 S data we assume optically thin emission to obtain a lower limit of the column density and, as a result, also for the abundance. 



TABLE 9 

Pipe Nebula core general properties with respect to core 109. 





Diameter 


Mass 


Nr 2 


«H 2 


P% a 


SPA" 


X(N 2 H + ) 


X(CN) 


X(C 3 H 2 ) 


X(CS) 




(pc) 


(M ) 


(10 2l cirT 2 ) 


(10 4 cirT 3 ) 


(%) 


(°) 


(10-") 


(io- u ) 


(it)- 11 ) 


do-' 1 ) 


Core 109 


0.063 


4.00 


47.60 


36.57 


11.0 


3.9 


2.20 


8.54 


20.6 


40.6 


Relative values 


Diameter 


Mass 




"Hi 


p% 


<5PA 


X(N 2 H + ) 


X(CN) 


X(C 3 H 2 ) 


X(CS) 


Core 109 


10.0 


10.0 


10.0 


10.0 


10.0 


10.0 


10.0 


10.0 


10.0 


10.0 


Core 40 


16.5 


6.3 


2.3 


1.4 


4.2 


22.2 


20.8 


49.1 


7.5 


16.4 


Core 14 


11.3 


3.5 


2.S 


2.5 


1.8 


40.4 


4.0 


11.2 


1.7 


68.2 


Core 48 


20.2 


5.2 


1.3 


0.6 


1.8 


83.8 


<2.5 


<2.6 


<0.6 


53.9 



(a) Franco et al. 12010) 

and the more intense emission, in particular of late-time 
molecules. Interestingly, the CCS abundance in core 109 is 
probably lower than in core 40, which is consistent with the 
view that the CCS molecule is destroyed soon after the forma- 
tion of a dense core, probably as a result of the contraction of 
the core dde Gregorio-Monsalvo et alj|2006t iMillar & Herbsll 
[T990tlSuzukietal.lll992h . This would reinforce the view that 
this core is in a very advanced evolutionary state. 

In summary, Core 109 seems to be the more chemically 
evolved core, probably because it is more dense and because 
it shows higher abundances of late-time molecules. Core 40, 
with three times lower column density, also shows large N2H + 
abundances. It might be in an intermediate chemical evolu- 
tionary stage. These two cores probably are in an evolution- 
ary st age slightly younger than that of the pr ototypical starless 
cores dTafalla et al.l2004blCrapsi et al.l2005h . Cores 48 and 14 
show similar physical properties in terms of size, mass and H2 
column density, to Cores 109 and 40. However, they appear 
to be very chemically poor and, therefore, they could be in an 
even younger stage of chemical evolution. 

4.4. Evolutionary trend and correlation with the diffuse gas 

Table [9] shows the summary of the main properties of the 
cores relative to Core 109, which is the one that shows the 
strongest line emission. In this table we show the physical and 
chemical properties. Additionally, we added the averaged po- 
larimetric properties of the diffuse envelope around the cores 



dAlves et al J 12008b iFranco et al.ll2010T) : polarization fraction 
(p%) and dispersion of the polarization position angle (SPA). 

As shown in Fig. 12 the polarization vectors calculated from 
optical extinction cannot be derived at the more dense regions, 
where the visual extinction is higher. In Fig. 12 except for the 
map of Core 48 with the lowest RMS, the polarization vectors 
lie in regions below the 3-cr noise level. However, the trend of 
the polarization vectors is in general rather uniform over the 
whole map. Indeed, there are vectors up to very close to the 
dense parts of the cores. Consequently, the derived magnetic 
field properties of the diffuse surrounding medium are also 
representative of those of the dense part of the cores. 

A relationship between the magnetic and the chemical 
properties of each core seems to exist. The two more chem- 
ically evolved cores, 109 and 40, appear to be embedded in 
a strongly magnetized environment, as SPA values clearly 
reflect (see Table |9}. The other two cores, 14 and 48, do 
not show very different morphological properties with respect 
to the previous two (size and mass). However, their chem- 
ical properties are completely opposed, and they are likely 
younger cores in chemical time-scale. Interestingly, the mag- 
netic properties of Cores 14 and 48 are also opposed to those 
of Cores 40 and 109. Cores 14 and 48 are surrounded by a 
molecular diffuse medium that is much more turbulent than 
that surrounding the two previous ones. Core 14 is possibly 
affected by the star formation undergoing in the nearby region 
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B59. Core 48 appears to be dominated by turbulence and con- 
stitutes an exception in the stem, whose cores have uniform 
mag netic properties amo ng them, showing low p% and high 
SPA dFranco et al.ll2010l) . 

In summary, these four cores of the Pipe Nebula have 
similar masses and sizes, but they are in different stages 
of chemical evolution: Cores 109 and 40 are much more 
evolved chemically than Cores 48 and 14. The different mag- 
netic properties of the diffuse molecular environment suggest 
that Cores 109 and 40 have grown in a more quiescent and 
slowly way (probably through ambipolar diffusion), whereas 
the growth of Cores 14 and 48 has occurred much faster, an 
indication that possibly a compress ion wave that generates 
turbulence or the turbulence i tself dFalle & Hartquistl 120021 : 
iBallesteros-Paredes et al.l2007l) . The longer time-scale of the 
ambipolar diffusion process could explain the more evolved 
chemistry found toward the cores surrounded by a magnetized 
medium. These features suggest two different formation sce- 
narios depending on the balance between turbulent and mag- 
netic energy in the surrounding environment. The importance 
of these results is worth of a more detailed study of the Pipe 
Nebula cores in order to fully confirm these trends. 

5. SUMMARY AND CONCLUSIONS 

We carried out observations of continuum and line emis- 
sion toward four starless cores of the Pipe Nebula spread out 
along the whole cloud selected in base of the ir magnetic prop- 
erties dAlves et a l. 2008; Franco et al.ll2.010t) . We studied their 
physical and chemical properties, and the correlation with the 
magnetic field properties of the surrounding diffuse gas. 

1 . The dust continuum emission of the observed Pipe Neb- 
ula cores shows quite different morphologies. In the 
sample there are diffuse cores, such as Cores 40 and 48, 
and compact and dense cores, such as Core 109. We 
have also mapped a clumpy filament, which contains 
the embedded Core 14. This filament is possibly under- 
going fragmentation into smaller cores of sizes com- 
parable to that of the others. We derived average radii 
of -0.09 pc (-18600 AU), densities of ~1.3xl0 5 cirT 3 , 
and core masses of ~2.5 M . 

2. The dust continuum peak coincides within the errors 
with Ay peak derived from the 2MASS catalog. The 
continuum emission is more sensitive toward the dense 
regions, up to > 10 magnitudes for the densest cores. On 
the other hand, the diffuse emission is better traced by 
the extinction maps. The masses are in average ~3.4 
times smaller. 

3. We have observed several early- and late-time lines of 
molecular emission toward the cores and derived their 
column densities and abundances. The starless cores of 
the Pipe Nebula are all very young, but they present dif- 
ferent chemical properties possibly related to a different 
evolutionary stage. However, there does not seem to 
be a clear correlation between the chemical evolution- 
ary stage of the cores and their position in the cloud. 
Cores 109 and 40 show late-time molecular emission 



and seem to be more chemically evolved. Core 109 
shows high abundances of late-time molecules and it 
seems to be the more chemically evolved. Core 40 
has three times lower H2 column density than that of 
Core 109. It presents a large N2H + abundance and the 
largest CN abundance, thus it might be in an interme- 
diate chemical evolutionary stage. Cores 48 and 14 
show only early-time molecular emission, and Core 14 
presents weak N^HT emission, and seem to be chem- 
ically younger than the other two cores. Core 14 has 
a similar mass and size than Core 40, but the N2H + , 
C3H2, and CS abundances are about one order of mag- 
nitude lower than the Core 40 abundances. Our re- 
sults and interpretation of the evolutionary stage of each 
core are consistent wi th the previous observations of 
iRathborne et al.l (120081) in these same cores. 

4. There seems to be a relationship between the properties 
of the magnetic field in the cloud medium of the cores 
and the chemical evolutionary stage of the cores them- 
selves. The two more chemically evolved cores, 109 
and 40, appear to be embedded in a strongly magne- 
tized environment, with a turbulent to magnetic energy 
ratio of 0.05 and 0.27, respectively. The two chemi- 
cally younger cores, 14 and 48, appear to be embed- 
ded in a more turbulent medium. This suggests that 
the magnetized cores probably grow in a more quies- 
cent way, probably through ambipolar diffusion, in a 
time-scale large enough to develop the richer chemistry 
found. On the other hand, the less magnetized cores 
likely grow much faster, probably in a turbulence dom- 
inated process, in a time-scale too short to develop late- 
time chemistry. 

5. The Pipe Nebula has revealed to be an excellent labo- 
ratory for the study of the very early stages of the star 
formation. The studied cores show different morpholo- 
gies, chemical evolutionary stages and magnetic prop- 
erties. The physical and chemical properties are not di- 
rectly linked as the competition between the magnetic 
field and turbulence at small scales seems to have an 
important influence in the core evolution. The impor- 
tance of these results require a more detailed study of 
the chemistry and magnetic field properties of the cores 
to fully confirm these results. 
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APPENDIX 

A. CALCULATION OF COLUMN DENSITY AND MASS OF DUST EMISSION 
Radiative transfer equation and Planck function 
The intensity emitted by an assumed homogeneous medium of temperature T ex and optical depth r v at frequency v is given by 

I v = B v (T et ){l-e^), (Al) 

where B v is the Planck function, 

2hv 3 1 

B v = — — — — -, (A2) 

c 2 e hv/kT _ ^ 

and c is the speed of light, k the Boltzmann's constant and h the Planck's constant. The Rayleigh- Jeans limit, hv «: kT (in 
practical units [v/GHz] « 20.8 [T/K]), does not hold for MAMBO-II observations (250GHz) of prestellar cores (r * 10K), 
preventing the use of this limit simplification. 

Telescope measurements 

The beam solid angle is = J^ eam PdQ. , where P is the normalized power pattern of the telescope. Assuming that the 

telescope has a gaussian beam profile, P reads P{6) = exp(-41n20 2 /6 , ^ pBW ), where is the angular distance from the beam 
center. The beam solid angle is 

° A = 4in72) ™- (A3) 
For discrete sources we measure flux densities, S v , instead of intensities, I v . These two quantities are related by 



S v = f 

This integration for a beam area, Sjj eam , allows us to calculate the beam-averaged intensity as 



IyPdQ.. (A4) 

source 



g beam 

<4) = -p— ■ (A5) 

From Flux to Column Density and Mass 

One can calculate the opacity of the emission measured inside a beam, rjf 11 " 1 , from equations I A 1 1 and I A5 1 relating it with the 
measured flux by 

/ o beam \ 

T beam = _ ln h _ — — L— — - . (A6) 



Qa B v {T) 



On the other hand, the optical depth is defined as 



= I K v pds, (A7) 

«7line of sight 



where k v is the absorption coefficient per unit density. 
One can relate the column density with the optical depth and, thus, with the measured flux using 



Nh 2 = I «h, ds = I -?—ds = Ty , (A8) 
p inn p inn k v 



which particularized to a beam is A^ 2 ™ = ^ naK „ where niu is the hydrogen mass and p is the mean molecular mass per hydrogen 
atom. In the case of optically thin emission the intensity is proportional to the column density as Eq. IA1I can be simplified to 
l v « B v (T)t v . 
Then, the mass can be calculated as 



l D 2 j N H2 d£l, 



M = pm H D 2 | N Hl d€l, (A9) 

which for a beam is M beam = p mu D 2 N^ m Qa, where D is the distance to the source. All these calculations can be applied to 
any solid angle bigger than a beam. 
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B. CALCULATION OF COLUMN DENSITY OF LINE EMISSION 

The column density for a / — > J - 1 transition of a molecule ('Mol') is 



NmoI - 



3 k Q mt e 



Eu/T a 



1 



Jy ( ^ex ) 



8^ 3 gK gi v s j U n 2 j v (r ex ) - j v (r bg ) 



1 - e T Jli„e 



T M B dv, 



which translates into useful units as 



(Bl) 



Mol 



cm 



= l. 6 7 X 10 14 -^i 



S 



Jkl 



gK gi [ erg cm 3 statC cm 2 



-2 



„Eu/T a 



GHz 



J v (T ex )-Jy(T bg )l-e- T 



Kkms 



(B2) 



Here, J v is the energy in units of temperature, and it reads J v (T) = e t% k _ x ■ See next subsections for a detailed description of all 
the terms involved. 

Observational terms 
Single transitions 

In case of single transitions, we have performed a Gaussian fit to the spectrum or a statistical moment calculation, both using 
tasks from the CLASS package. We obtain from either analysis the main beam temperature, Tmb, the line velocity, v, and the 
integrated emission, J Tmb dv. 

The opacity, t, is calculated numerically in those molecules with more than one transition observed. In the other cases we have 
assumed t ~ 0.3. The excitation temperature, r ex , can be calculated from the radiative transfer equation as 



hv 

T - 

1 ex , 

k 



In 



hv/k 



+ 1 



(B3) 



where 7b g is the background temperature. 



Hyperfine transitions 

In case of hyperfine transitions, we take into account all the hyperfine components of the selected transition. We have performed 
a hyperfine fit using CLASS, which provides A x t,„ , v ^ encellne , Av , t,„ , where A is 



A = f(j v (T eK )-J Y (T hg j), 



(B4) 



beeing / is the filling factor assumed to be ~1. 



To be able to use equation lB2l as in the single transition case, we need r ex , t and J Tmb dv. We can calculate T ex as in equation 
"calculating Tmb as A x T m /r m , and t,„ is given by CLASS. For the integrated emission, we can use 



to Av 



-f 

«71ine 



rdv ^ 



1 



To 



J v (T ex ) - J v (Tb g ) 1 - e- 



f 

«71im 



?mb(v) dv, 



leading to 



J 7mb(v) dv a t Av (j v (T ex ) - 7 v (r bg )) 

J line 



1 - e- 70 



TO 



(B5) 



(B6) 



Making this transformations eq. IB2l can be used for hyperfine transitions. 

Non-observational terms 
Partition function (Q YOt ) 
The rotational partition function, Q m t{T), is defined as 

QUT)^JS k gie- hBJ(]+i)lk \ 



(B7) 



where the gx factors are the degeneration of the respective quantic number, in particular gj = 2 J + 1 . 

Eq. IB7l can be approximated, in the limit of high temperatures, by an integral because generally the energy levels are close 
together. We are only interested in the high temperature limit because is when the transition is activated, so this limit is accurate 
enough. 
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• Linear molecules: The solution for the diatomic case is general for any lineal molecule, so long as the molecular moment of 
inertia is computed properly for more than 2 atoms. 

For lineal molecules g^ = 1, gi = 1 and gj - (27 + l)/cr. <x (the symmetry number) is 1 for heteronuclear diatomic (C-O) 
or asymmetric linear polyatomic (O-N-N) molecules, and 2 for homonuclear diatomic (H-H) or symmetric linear polyatomic 
(O-C-O) molecules. 

The partition function at high temperatures can be calculated as 



Qn 



i r 

o~ Jo 



(2/+l)e 



-hBJ(J+l)/kT 



dJ 



Q -(J 2 +J)hB/kT 



d(J 2 + J)) 



1 kT 
v~hB' 



(B8) 



where B is the rotational constant available at the catalogues. A more accurate expression (Pickett et al., 2002) used in this work 
is 



1 (kT 1 



Qxot *a\hB + 3 



1 crhB 
15 kT 



+ 



(B9) 



• Non-linear molecules: Non-linear molecules have up to three moments of inertia and, thus, three rotational constants (A, B, C). 
In a similar way than before, but more complicated, the calculation of the rotational partition function at high temperatures is 



Qr, 



£(¥) 



3/2 



1 



VaTJc 



(BIO) 



Upper level energy (E„) 

We can calculate the energy of the upper level (E u ) as a function of the lower level (£/) plus the energy of the photon emitted 
(both available at catalogues). This is, in units of temperature and using the units given in the catalogues, 



Eu 
K 



1.4388 



+ 4.799 x 10" 



MHz 



(Bll) 



Intrinsic line strength times squared dipolar momentum (S fi 2 ) 
We can calculate the product of the intrinsic line strength, S jui, and the squared dipolar momentum, /i 2 , from the Q mt at 300 K 



(Grot )' t ' le l me strength (LoglNT) at 300 K and the lower state energy (£/). All these parameters are available at the catalogues. 
In a usable form, 



erg cm 3 statC 2 cm 2 D 



\-2 



24025 x 



\QLogINT 



MHz nm 2 



Q 



300 



MHz 



exp 4.796 x 10" 



Ei 



cm 



)(l-exp{- 



1.6 x 10" 



MHz J 



(B12) 



